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In a first order approximation, the Bose-Hubbard Hamiltonian with on- 
site interaction is obtained from the free Hamiltonian [U = 0) and generalized 
commutation relations for the annihilation-creation operators. Similar gen- 
eralized commutation relations were used for the first time in high energy 
physics. The spectrum of the system can be found formally by using the 
algebraic properties of the generalized operators. 
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In recent years much attention has been given to the Hubbard model for 
bosons ([1-5] and references therein). At first glance, this is a simple descrip- 
tions for interacting particles, and related to a subject like phase transitions. 
In this brief report we obtain, in a first order, the Bose-Hubbard model from 
usual phonon-like Hamiltonian where the annihilation and creation opera- 
tors satisfy generalized canonical commutation relations. In fact, the Bose- 
Hubbard model corresponds to the limit of small perturbation with respect 
to usual commutation rules. To the best of our knowledge this connection 
is not found in the literature. Curiously, the model is formally solvable in 
a way similar to the frcc-phonon system. So, this work opens new scopes 
about the interpretation of the Bose-Hubbard Hamiltonian, and its relation- 
ship with some ideas in high energy physics. The feeling of the report is that 
interaction terms can be considered modifying the commutation rules for the 
relevant operators. 



Consider the Hamiltonian H, and the operator M defined by (we use units 
in which U = 1) 

n^jY^MPlPk + PkPl), (1) 
^ k 

^f-^J:iPiPk+pkPl), (2) 

k 

where we assume ujk a real function of the wavenumber k (first Brillouin 
zone) and L is the total number of sites in the lattice. Moreover, consider 
the operators /3k and /3l satisfying the implicit commutation relations 

[P,,Pl] = (L + -M)6k,,. (3) 

^k 

where 5k,p is the usual Kronecker symbol and C/ is a constant which will 
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be related to the usual strength-interaction-parameter in the Bose-Hubbard 
model. Moreover we assume 

[f3k,Pp]^[PlPi]^0. (4) 
With respect to the above definitions we remark : (a) Similar general- 
izations of the canonical commutation relations were proposed originally in 
high energy physics to explain some phenomena related to quark physics [6- 
8] where the correction was proportional to the Hamiltonian rather than to 
the operator Af . (b) The case [7 = corresponds to the well-known model 
of free excitations with spectrum Uk, where PI and Pk are the usual bosonic 
operators, (c) We can show [Ti^Af] — [9] and then, the commutation rela- 
tions (3) do not change with time, (d) If we consider the limit L — > oo then 
Y J2k ~^ J dk in (1) and (2). Moreover, in this hmit we have L5k,p — > S{k—p) 
and a correction proportional to the density ^^e commutator (3). 



To determine the first order relationship between the system defined by 
(1-3) and the Bose-Hubbard model we go on as follows: let bk and bl be the 
usual creation- annihilation Bose operators which commutator [bk, &p] = LSk^p. 
Since U ^ we expect (3k —>■ bk ; then, it seems reasonable to assume the 
first-order expansion in U, 

h^bk + ^Y. f(k, P, r, s%b% + 0{U^) (5) 

^^fc prs 

where / is a unknown real function which will be determined using, at first 
order, the generalized commutation rules (3) (i.e. we put (5) in (3)). After 
some calculations we get 

Pk^hk + E k^r,p^sbpb% + 0(C/2) (6) 
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Namely, the expansion (5) is consistent with the commutation rules (3). An 
homogenous zero order groundstate was assumed. In this approximation (6), 
the Hamiltonian (1) becomes 

n = \T.^Mh + hbl) + ^^ ^k+r,,+s (hblbrbl + bib,b%)+0{U'). (7) 
fc krps 

And if we consider the Fourier transform to real space / in the lattice 

k^y^bke^'^ (8) 
^ k 

then, the Hamiltonian (7) can be written at first order in U as 

^ = - E t(blbi_, + bibl, + bibU + blk+,) + U J2(blbMk + bibjbibl) + 0{U^) 
I I 

(9) 

which corresponds to the symmetrized Bose-Hubbard Hamiltonian of strength 
U in position representation. To obtain (9) the usual dispersion relation in 
one-dimension uj^ — —2tcosk, where t is the hopping parameter, was as- 
sumed. It must be noted, however, that the Hubbard term is also obtained 
for dimension greater than one {D > 1). 

To determine the spectrum of the system (1-3) we consider the eigenvector 
of n and J\f given by \Ek,n; Nk,n)- Namely, 

I Ek^n', Nk,n) — Ek,n \ Ek,n', A^it,n) , (10) 
•A/" I Ek^n]Nk,n) = Nk,n \ E^^n, Nk,n) • (H) 
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To find the relationship between different eigenvalues of Ti (i.e. between 
Ek,n+i and Ek,n) we go on as follows : by using the commutation rules (3), 
one finds 

n(3l\Ek,n, Nk,n) = {Ek,n + ^LOk + T + Pk^)}\Ek,n\ Nk,n) ■ (12) 



However, also from equation (3) one has 



1-7^ 

Lu>k 



(13) 



then f^l I Ek^n'i Nk^n) is an eigenvector of J\f. Expressions (13) gives the 
relationship between different eigenvalues Nk^n of the operator J\f : 



2 + ^M(l + i£r) 



k,n+l 



1-7^ 

Luik 



(14) 



Introducing (13) onto (12), we have Hf^k \ ^fc,„; Nk,n) = Ek,n+i/3k \ Ek^n, Nk,n 
and the relationship for the spectrum of Ti. 

2U 



Ek,n+i — Ek,n + 2a;fc + 



L- 



u_ 



iNk,n + l) 



(15) 



where if [/ = becomes the usual oscillator-likc-spcctrum for phonons. The 
spectrum (15) is not given in explicit form because of its dependence on the 
eigenvalue Nk,n ■ Nevertheless, from equation (14), it can be solved if we 
note that 



Nk,n^ 



Luk 



U 



1 



u 

Luju 



Luik 



(16) 



where it was assumed = 0. 
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We finally raise some questions related to the above connection between 
the Bose-Hubbard Hamiltonian and the generalized commutation rules: (i) 
The extension to the usual Fermi-Hubbard model was not worked-out. The 
generalization of the commutation (or anticommutation) relations (3) is not 
straightforward, (ii) Thermodynamic properties like internal energy, heat ca- 
pacity , thermal conduction, must be calculated by using the spectrum (15). 
Nevertheless, this seems not direct, although one may attempt a perturba- 
tion expansion in powers of U. (iii) The usual term related to density-density 
nearest-neighbor interaction [2-5], between bosons, is absent in (9). It seems 
that such a generalization requires a modification of the commutation rela- 
tions (3). (iv) In last years, much attention has been given to disordered 
Hubbard models ([1, 10-13] and references therein), in our case it is not clear 
how to consider disorder in the formulation (1-3). Nevertheless, this can be 
partially carried-out using an appropriate spectral function Uk like to this 
one of disordered systems. 

A more detailed treatment of points (i-iv) and further physical applica- 
tions will be given elsewhere. 
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